Incommensurate spin correlations in Heisenberg spin-1/2 zig-zag ladders 
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We develop a low-energy effective theory for spin-1/2 frus- 
trated two-leg Heisenberg spin ladders. We obtain a new type 
of interchain coupling that breaks parity symmetry. In the 
presence of an XXZ-type anisotropy, this interaction gives 
rise to a novel ground state, characterized by incommensurate 
correlations. In the case of a single ladder, this state corre- 
sponds to a spin nematic phase. For a frustrated quasi-one- 
dimensional system of infinitely many weakly coupled chains, 
this state develops true three dimensional spiral order. We 
apply our theory to recent neutron scattering experiments on 
Cs 2 CuCl 4 . 

PACS No: 75.10.Jm, 75.40. Gb 



Quantum spin chains have for a long time attracted the 
attention of both theorists and experimentalists. One of 
the main reasons for this continuing fascination is the 
dominant role played by quantum fluctuations in these 
systems, which lead to a rich variety of observed physical 
phenomena. In recent neutron scattering experiments |l]] 
on the quasi one-dimensional frustrated Heisenberg anti- 
ferromagnet CS2CUCI4, an intriguing type of spiral order 
was observed. The very existence of such spiral order, 
which is an incommensurate structure, is puzzling from 
a theoretical point of view. According to the standard 
lore, three dimensional ordering of quasi-lD systems re- 
sults from stabilization of the dominant spin correlations 
in the underlying ID constituents. Therefore, 3D spiral 
order would require strong incommensurate ID spin cor- 
relations. This is in contradiction with the known prop- 
erties of simple antiferromagnetic Heisenberg chains and 
ladders, where the only known mechanism for generating 
incommensurabilities is via external magnetic fields. 

In the present work we propose a novel mechanism 
that naturally gives rise to incommensurate correlations 
in spin ladders and quasi ID materials in the absence 
of external fields. Interestingly, we find that this phe- 
nomenon occurs in a standard model of two coupled spin- 
1/2 zig-zag Heisenberg chains, which has attracted much 
recent interest H,||. 

In the bulk of this letter we shall concentrate on this 
simple frustrated ladder system in order to clearly exhibit 
the precise nature of the mechanism of incommensurate 
spin correlations. The application of our results to the 
quasi ID case relevant for CS2CUCI4 is briefly discussed 
at the end. In what follows we derive the low-energy 



effective theory for the zig-zag ladder system, which we 
find to contain a parity-breaking interchain interaction. 
This "twist term" , which was not considered in previous 
studies, is an operator of conformal spin 1 Q and has 
important consequences: We will show that, although it 
alters the position of phase boundaries already in the spin 
rotationally invariant (SU(2)) case, it leads to incommen- 
surate correlations only in the presence of an easy plane 
XXZ anisotropy. 

The Hamiltonian of the anisotropic zig-zag Heisenberg 
ladder is 

H = J X! y"X S J> ' s j>+i + (A - l)Sj >n Sl n+1 } 

j=l,2 n 

+ J' J]] + Si !TS _|_i} • S 2 ,„ 
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where \J'\ -C J are the exchange couplings, and 
1, 1 A' I < 1 are the anisotropy parameters. 

In order to derive the low-energy effective field the- 
ory for (|l|) we use standard bosonization techniques [Q. 
In the continuum limit the spin operators decompose as 
Sj,n - * a o [Jj( x ) + {~l) nn j{ x )]i where ao is the lattice 
spacing and x — nag- Here 3j(x) and iij(x) are the 
smooth and staggered components of the magnetization 
operator. In the framework of the bosonization method 
these quantities are expressed in terms of canonical Bose 
fields $j(x) and their dual counterparts Qj(x), where 
[&j(t, x),$j(t, x')] = —i9(x — x 1 ). Changing to symmet- 
ric/antisymmetric combinations $± = ($i±<I>2)/v / 2, the 
Hamiltonian density consists of several terms. The "free" 
part is 



(2) 



<j=± 



where v a oc Jclq are the spin velocities. In addition, there 
is an in-chain current-current perturbation 



Hjj — 9i cos %/47r$ + cos V47r$- 
+52 V [{d x Q a ) 2 - {d x $ B 



<T = ± 



(3) 



The interchain interaction gives rise to two perturbations 
which we denote by 1i.ee and Tips respectively 



1 



TLcc — 53 cos V4vr < i > + cos \/47r0_ 

-, 94 ° [(^e ff ) 2 - (dM 2 } , (4) 

Kpb = ff5(-9a!$- sin\/47r$ + + d x § + sin\/47r$_) 

+56^6+ sin -\/47r0_ . (5) 

The perturbation TLcc is the well-known current-current 
interaction |^]||, which promotes dimerization and leads 
to the formation of a spectral gap. The perturbation 
TLpb is the novel parity breaking term, which in terms 
of the staggered magnetizations 711,2(2;) reads 

Hpb ~ ni • d x n 2 - n 2 ■ 9 x ni 

+ (A' - 1) \n\d x nl - nld x n\\ . (6) 

The in-chain coupling constants 51.2 are determined by 
J and A, whereas gs,...,ge are functions of J', A'. 

Around the SU(2) symmetric point A = A' = 
1, the bosonized Hamiltonian can be expressed in 
terms of four Majorana fermions [§,1). The Hamilto- 
nian of two decoupled chains simply becomes TLo = 

" l f E?=o (&Mr ~ • The perturbation H> = 

TLjj + TLcc +TLpb is of the form 

3 4 

ff' = J>i Aj+YthBj, (7) 

j=i i=i 

where the operators Aj and are given by 
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We note that in the above formulas we have neglected 
terms that lead to a renormalization of the spin velocities 
H . The couplings aj , (3k are easily expressed in terms of 
the gi. The operators Aj originate from TLpb, whereas 
the -Bfc's stem from TLcc and TLjj. 

All operators are marginal but while the have con- 
formal spin 0, the Aj have conformal spin 1. We de- 
rived the renormalization-group flow for the perturba- 
tion TLpert from the Operator Product Expansion for the 
perturbing operators in the spirit of |?|,[| 

Q'i = 2a 2 /3 2 + 4a 3 /3i 

d 2 = 2ai/3 2 + 4a 3 /3 4 

a 3 = 2ai/?i + 2a 2 /3 4 + 2a 3 /3 3 

/9i = -4a 2 a 3 + 2/3i/3 3 + 2/3 2 /3 4 (9) 

02 = -4ai + 4/3i/3 4 

$3 = -4a ia2 + 2f3j + 201 

fa = -4a ia3 + 2/3i/3 2 + 2/3 3 /3 4 



Here a dot denotes the derivative with respect to the RG 
logarithm const x ln(Au s /|a->|), where A is a momentum 
cutoff H . The RG flow determined by (||) is clearly com- 
plicated. Therefore we performed a numerical analysis of 
(0). Our findings are the following: 

• At the SU(2) symmetric point, the current-current 
interactions reach the strong coupling regime first (in the 
cases where the flow is towards strong coupling). This 
corresponds to the dimerized phase found in ||. A new 
feature caused by the twist operators is that the tendency 
towards dimerization extends into part of the region of 
ferromagnetic current-current interchain interactions (see 
Fig. @) 
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FIG. 1. RG flow in the SU(2) symmetric case. There 
are only 2 independent couplings for twist (01,2) and cur- 
rent-current operators (/3± = (03 ± /3i)/2) respectively. The 
initial conditions are chosen such that the couplings f3± (0) < 
and in the absence of the twist operators would flow to zero 
coupling. 

• Away from the SU(2) symmetric point both twist and 
current-current couplings generally (but not always) flow 
towards strong coupling. However, it is now possible for 
the twist terms to reach the strong coupling regime first 
(see Fig. ||)! We take this as evidence for the existence 
of a novel phase, the physics of which is determined by 
the twist terms. 

In order to elucidate the physics of this new phase we 
now turn to a mean-field analysis of the perturbations. 
The situation is particularly simple in the limit of very 
strong anisotropy, |A|,|A'| w 0, corresponding to two 
coupled XX chains. We therefore shall discuss this case 
first and only then generalize to arbitrary anisotropy. Us- 
ing the bosonization formulas for the XX point |||| an d 
retaining only the relevant part of the perturbing opera- 
tor we arrive at 

TL = TL +79 x e + sin\/27e_, (10) 

where TLo is given by (||). 
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FIG. 2. Example of an RG flow where the twist couplings 
reach strong coupling first. The structure of equations (^) is 
such that the twist couplings cannot grow without "help" of 
the current-current interaction. 

In order to analyze ( |To| ) we use a self-consistent mean- 
field approach. Assuming that the ground state of the 
system is found in the sector with a nonzero topological 
spin current d x Q+, we arrive at the following mean- field 
Hamiltonian 



Umf =Ho + Kd x Q + -/iAsinV27r8_, (11) 



where 



7 (sin\/2^e_), /iA = - 7 (9 x e+). (12) 



Thus Umf decomposes into 2 commuting parts, Ti+ 
Ti.-, with 



H- = y ) 2 + (^e_) 2 ] - ^Asin^e. 



H+ 



(&e H 



2k 



(13) 



The "+" channel is solved by eliminating the d x Q+ term 
through a field redefinition: 0+(x) — ► 0+(.t) — kx/v s - 
The average value of 9^0+ is then given by (d x ®+) = 
— k/v s . The "— " channel is a Sine-Gordon model for 
the dual field and can be solved exactly, the expec- 
tation value of the mass term being (sin y27r©— ) = 
c(\(j,\/ Avs) 1 ' 3 sgn fj,, where the constant c can be calcu- 
lated JnJ. The self-consistency conditions (]l^) then lead 
to the solution 



/i = ±Av s z c 



K = ±Av s z c 



2 3/2 



(14) 



where zrj = j/Av s is a dimensionless coupling con- 
stant. The expectation value of 0_ is determined by 
the position of the minima of the Sine-Gordon poten- 
tial; (©-) = ^/7r/8 sgn /x(mod \/2tt) . This allows us to 
express the dual fields in chains 1 and 2 as 



e li2 (x) = ^©°+(z) 



V2v t 



± ^ , (15) 



where @°,(x) as well as its derivative have zero expecta- 
tion value. Using this together with the bosonization ta- 
ble at the XX point, we find the following asymptotic be- 
haviour for spin-spin correlation functions in the strong- 
coupling phase 



(St(x)S-(O)) 



11/4 



exp 



^/2(kx/v s ) (16) 



where j = 1,2. The transverse spin correlations still fall 
off with a power law but their decay is much slower than 
for an isolated XX chain where (n + (x)n~ (0)) ~ |x| -1 / 2 . 
Furthermore the correlations are incommensurate] The 
characteristic momentum of the magnetic spiral is 



q Q = 7r/a - 



(17) 



-• L 



The deviation from the anitiferromagnetic wave vector 
7r/ao is very small for a weak interchain coupling, which 
is in qualitative agreement with Q. We have performed 
a stability analysis of the above mean-field solution and 
found only convergent corrections, e.g. to the exponent 
1/4 in 

The solution of the XX zig-zad ladder ([n]) we obtained 
actually describes a spin nematic ground state of the 
model. This phase, with unbroken time reversal sym- 
metry pdf| is characterized by nonzero local spin currents 
polarized along the anisotropy (z-)axis. The longitudinal 
(in-chain) component of the total spin current is given by 
Ju(x) = —y / 2/Trv s d x Q + (x). Using equations of motion 
for the spin densities, one easily finds the transverse (in- 
terchain) part of the current, originating from the twist 
term in @: Jf_{x) = 



From the 

mean-field analysis above one immediately obtains the 
important result that (J,f) = —{Jj) = ±^2/^7. Thus, 
in the ground state the longitudinal and transverse spin 
currents are equal in magnitude but propagate in oppo- 
site directions. The resulting picture shown in Fig. 3 
demonstrates local currents circulating around the trian- 
gular plaquettes in an alternating way, with the total spin 
current of the system being zero. 

The spin nematic phase preserves the spin U(l) sym- 
metry but spontaneously breaks a special discrete Zi 
symmetry of the model. Indeed, the twist term (|^) 
is invariant under a tensor product of the site-parity 
and link-parity transformations Jl2| on the two chains, 
V 12 = Pf ® Pi [though it breaks Pf (L) ® P 2 5(L) ]. In the 
XX case (|l0|), this transformation reverses the signs of 
all currents, reflecting the Z2 degeneracy of the ground 
state. The spectrum of the system is rich and, apart from 
gapless excitations in the " +" channel, contains massive 
quantum solitons and their bound states (breathers), as 
well as kinks carrying fractional topological charge. 
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FIG. 3. Structure of the spin currents in the spin nematic 
phase. 

Moving in parameter space from A, A' = (XX point) 
to A, A' = 1 (XXX point) in the general XXZ case, we 
found the following: 

• The spin-nematic phase occupies the whole region 
<?4 < 92 beyond which a gap appears in the + channel 
as well, and the spin correlations become short-ranged 
(although they remain incommensurate). 

• In the fully gapped phase, the spin currents and the 
incommensurate wave vector go & r e reduced and continu- 
ously vanish as the SU(2) symmetric point is approached. 
Indeed, the spin nematic state cannot be stable at the 
XXX point because nonzero spin currents would break 
the SU(2) spin rotational symmetry, which is forbidden 
in 1 + 1 dimensions. 

In light of the experimental results yj] it is important 
to investigate the effects of an external magnetic field. At 
the XX point we find that a longitudinal field has essen- 
tially only trivial effects, whereas a transverse field leads 
to the formation of a gap in the + channel and completely 
destroys the spiral phase (above a critical field). Indeed, 
the magnetic field oriented along the x-axis in spin space 
enters the Hamiltonian with the term |l3fl 

= h ^2 sin V47r$j cos y/nQj 
i=i,2 

where h is proportional to the applied field. This pertur- 
bation again has conformal spin 1 and, together with 
the twist term (jU^), leads to a complicated RG flow, 
where an effective XY anisotropy, described by the fol- 
lowing conformal spin zero operator, is generated p3| : 
Oxy ~ h 2 cos \/2ttQ+ cos v2tt©-- This strongly rele- 
vant operator locks the fields 0± at the vacuum values 
0+ = \Ar/2, 8_ = 0, or vice versa, and gives rise to a 
commensurate Neel ordering of the spins along the y-axis, 
characterized by a nonzero average staggered magnetiza- 
tion (n\ + n\) ~ (sin ^72e + )(cos ^726-). 

Let us now turn to the case relevant for the exper- 
iments where we have two-dimensional arrays of 
weakly coupled chains. We illustrate our findings for the 
simpler XX case. The perturbation to the free part of 
the multi-chain Hamiltonian is 

3 



We perform a mean-field analysis in a way similar to 
the two-chain case above. Assuming that the following 
averages 7(sin y/n(Qj - Qj+i)) = k, 7<9 x (@j + ©j+i)) = 
2\x are nonzero and redefining the fields as Qj — —2nx + 
jy/ir/2 + Qj , with (Qj) — we obtain a self-consistent 
mean-field Hamiltonian 

Hmf = n [Qj] - 2(mJ2 C0S V5F(©j - Qj+i) ■ (19) 

3 

The Hamiltonian ([l9]) can be viewed as describing cou- 
pled Josephson-junction arrays and leads to the pinning 
of the fields Qj . The resulting average staggered magne- 
tization is given by 

(n ± (a;)) = exp(±wrj/2 ± 2^nx) . (20) 

This corresponds to incommensurate spiral order along 
the chains with a 90° rotation of the average staggered 
magnetization in the transverse direction. This is in qual- 
itative agreement with experiment Q. Details of the 
above calculations as well as a more quantitative com- 
parison to experiment will be presented in a separate 
publication j|. 

We thank P. Azaria, R. Coldea, R. Cowley, T. Gia- 
marchi, P. Lecheminant and A.M. Tsvelik for important 
discussions. 
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